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Abstract: A theoretical study of differences in broadband high-index- 

contrast grating (HCG) reflectors for TM and TE polarizations is presented, 
covering various grating parameters and properties of HCGs. It is shown 
that the HCG reflectors for TM polarization (TM HCG reflectors) have 
much thicker grating thicknesses and smaller grating periods than the TE 
HCG reflectors. This difference is found to originate from the different 
boundary conditions met for the electric held of each polarization. Due 
to this difference, the TM HCG reflectors have much shorter evanescent 
extension of HCG modes into low-refractive-index media surrounding the 
HCG. This enables to achieve a very short effective cavity length for VC- 
SELs, which is essential for ultrahigh speed VCSELs and MEMS-tunable 
VCSELs. The obtained understandings on polarization dependences will be 
able to serve as important design guidelines for various HCG-based devices. 
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1. Introduction 

The grating structure referred to as high contrast grating (HCG) is near subwavelength grat¬ 
ings formed in a thin high-refractive-index layer surrounded by lower refractive index mate¬ 
rials m. HCGs can provide high reflectivity to a surface-normal incident wave over a broad 
wavelength range, e.g. higher than 99.9% over broader than 100 nm. This capability makes 
HCGs an attractive alternative to conventional distributed Bragg reflectors (DBRs) for vertical- 
cavity surface-emitting lasers (VCSELs). Since the first demonstration of VCSELs employing 
a HCG reflector Eld, many novel HCG-based VCSEL structures have been reported, fea¬ 
tured by strong single-transverse-mode operation EHSl, efficient MEMS tuning of emission 
wavelength EUl, integration onto a silicon-on-insulator (SOI) wafer ISHD, and control of 
transverse mode confinement Gam. 

These HCG reflectors can be classified into TM HCG and TE HCG reflectors, depending on 
the polarization of an incident wave to which a HCG is highly reflective. Since both types can 
provide broadband high reflectivity and their reflection processes are explained by the same 
theories iniiBj, less attention has been paid to the differences between the two types. I have 
found from various TM and TE HCG cases ca that optimal HCG parameters such as grating 
thickness, period, and duty cycle, and the evanescent tail length of HCG modes into the low 
index media are quite different for each type. We need to note that this evanescent tail length 
can be a key parameter in the designs of HCG-based high-speed VCSELs and MEMS-tunable 
VCSELs, significantly influencing their performances as discussed in Section 4. This rather 
empirical observation and the applicational importance of polarization dependent properties 
have motivated the studies of this paper, asking two questions: whether the empirically observed 
polarization dependences occur by chance or result from a physical origin and whether other 
properties of HCG reflectors such as penetration depth, reflection phase shift, and stopband 
width are also polarization dependent. 

This paper is organized as follows: Eirstly, the polarization dependences of HCG parameters 
as well as HCG reflector properties are statistically investigated in Section 2. Then, the origin 
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Fig. 1. Field profiles around (a, b) TM and (c, d) TE HCG refiectors with 1550-nm- 
wavelength TM and TE waves being incident from the bottom, respectively. In (b, d), the 
field profiles at v=0 nm (tot) are shown together with their 0th (Oh) and 1st harmonic (Ih) 
components. The TM HCG has a grating thickness of 420 nm, a grating period of 642 nm, 
and a grating bar width of 410 nm while the TE HCG has a thickness of 245 nm, a period 
of 1084 nm, and a width of 322 nm. The black dashed lines denote grating boundaries. The 
refractive index of the grating made of Si is 3.47 and both incident and exit media are air. 


of polarization dependences are investigated in Section 3. In Section 4, the importance of po¬ 
larization dependent properties in designing HCG-based devices, mostly VCSELs is discussed. 

For numerical investigations, the rigorous coupled wave analysis (RCWA) method was em¬ 
ployed ini-i chose the design wavelength, grating material, and surrounding materials to be 
1550 nm, Si with a refractive index of 3.47, and air, respectively. However, the results in this 
paper should be valid for other wavelengths and other sets of materials with a high refractive- 
index contrast. 

2. Statistical investigation 

2.7. Insight from examples 

Figure 1 compares exemplary TM and TE HCG reflectors, both optimized to have a broad 
stopband around a wavelength of 1550 nm. In Figs. 1(a) and 1(c), the TM HCG reflector has a 
thicker thickness, a shorter period, and a larger duty cycle (ratio of grating bar width over grat¬ 
ing period) than the TE HCG. Regarding the evanescent tail, the TE HCG reflector has a much 
longer fleld extension into the transmitted side (above the HCG) than the TM HCG reflector, 
as shown in Figs. 1(a) and 1(c). Figures 1(b) and 1(d) decomposing flelds into harmonics show 
that the first order harmonic component Ih mainly determines the evanescent tail length and it 
extends not only into the transmitted side but also into the reflected side. 

2.2. Statistical analysis 

In order to And out whether the polarization dependences observed for the exemplary HCG 
reflectors are universal, a statistical investigation is conducted. A broad range of TM and TE 
HCG structures are simulated by varying the HCG parameters, i.e., grating thickness, grating 
period, and grating bar width independently over very broad ranges. Among the simulated 
parameter sets, only those HCG parameter sets that result in a reflectivity higher than 99.9 
% throughout a stopband broader than 100 nm around 1550-nm wavelength are counted for 
analysis. The histograms in Fig. 2 summarize the simulation results: 
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Fig. 2. (a)-(c) Histograms of HCG parameters (grating thickness, grating period, and grat¬ 
ing bar width), (d)-(f) Histograms of HCG properties (reflection phase, stopband width, 
and penetration depth). 
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Fig. 3. Plots of HCG properties (stopband width, reflection phase, and penetration) as a 
function of HCG parameters (grating period and grating width) at a speciflc grating thick¬ 
ness. The thicknesses of TE HCGs and TM HCGs are 245 nm and 420 nm, respectively. 


• Figures 2(a) to 2(c) show the polarization dependence of HCG parameters. As in the ex¬ 
emplary HCGs, the TE HCGs have thinner thicknesses and longer periods than the TM 
HCGs. The grating thickness and period values are concentrated around specific values 
while the the grating bar width values are widely distributed. Thus, one may guess that 
there could be conditions to be met for grating thickness and period, which are investi¬ 
gated in Section 3. 

• Figures 2(d) to 2(f) show the polarization dependence of HCG properties, i.e., refiection 
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Fig. 4. Transmittance plots of (a) TM HCGs and (b) TE HCGs as a function of wavelength 
and grating thickness. The grating period and width are the same as in Fig. 1. 


phase, stopband width, and penetration length (HKIHI. The TM HCGs are more likely to 
have a smaller reflection phase and a smaller penetration. Especially, the smaller penetra¬ 
tion leads to a shorter effective cavity length advantageous for high-speed VCSELs and 
MEMS-tunable VCSELs, as discussed in Section 4. Regarding the stopband width, TM 
and TE HCGs have similar distributions. 

To And out whether there is an apparent relationship between the HCG properties and HCG 
parameters, the HCG properties are plotted in Eig. 3 as a function of grating period and grating 
bar width at a given grating thickness: 

• As shown in Eig. 3(f), the penetration depth monotonically increases with the grating bar 
width. It can be explained in this way. As the grating bar width increases, the coupling 
between neighboring grating bars increases. Since the penetration depth is related to the 
lateral propagation distance within a HCG layer, the increased coupling may lead to a 
longer penetration length. The reflection phase also gradually increases with the grating 
bar width as shown in Eig. 3(e). 

• The bandwidth is hardly related to either the grating period nor bar width, as shown in 
Eigs. 3(a) and 3(d). 

• Other correlations shown in Eigs. 3(b) and 3(c) are not clear and sometimes change at 
different grating thicknesses. 


3. Origin of polarization dependence 

In this section, the origin of polarization dependences observed in grating thickness, grating 
period, and evanescent tail length are investigated. 

3.1. Grating thickness 

In Eig. 4, the transmittances of the exemplary TM and TE HCGs of Eig. 1 are plotted as a func¬ 
tion of wavelength and grating thickness. Broad bandwidths can be obtained at the bendings of 
low transmittance line ca. Eor both TM and TE HCGs, the convex bending occurring at the 
smallest thickness, designated by a green arrow gives the flattest curvature corresponding to the 
broadest bandwidth: other bendings with a thicker thickness do not provide a stopband broader 
than 100 nm. This explains why there is a single (not multiple) optimal grating thickness for 
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Fig. 5. Eigenmode electric-field profiles Et^m of (a, b) the TM waveguide array and (c, d) the 
TE waveguide array, with propagation constant, /3 values. The waveguide period and width 
are the same as given in the caption of Fig. 1. The black dashed lines designate waveguide 
boundaries. 


each polarization in Fig. 2(a). The reason this bending occurs at a shorter thickness in the TE 
HCG than the TM HCG can be explained by investigating the eigenmode profiles. 

Figure 5 shows the field profiles of the two propagating eigenmodes of the TM and TE 
waveguide arrays (along y axis) at a wavelength of 1550 nm, which are infinite versions (along 
z direction) of the TM and TE HCGs. There are only two propagating modes for both TM and 
TE cases as the considered wavelength is within the so-called two-mode regime ina . To have a 
transmittance close to 0, the lateral average of these two modes at the grating exit plane should 
cancel each other. This condition can be expressed by llT5]| . 

y (am + a^)A~^ [ Ei^m{x)e~‘^"’^’‘dx = 0, (1) 

m=l,2 Jo 

where is the transmittance amplitude of the 0-th harmonic component, am and are the 
coefficients of the upward- and downward-moving m-th modes at the exit plane, respectively, 
A is grating period, Et^m is the normalized electric field component for TM and Ey^m for 
TE) of the m-th mode, and Zex is the z coordinate of the exit plane. This Zex value becomes 
same as the grating thickness, if the input plane of the grating is assumed to be at z=0. Since 
the two propagating modes have similar amplitudes, i.e., a\ ^ a 2 and a^ a^, Eq. o can be 
simplified: 

rA 
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= {ai+aP)A-^ £ (2) 

m=l,2 

= (ai ^ 

where the normalization condition is used in the derivation from the first line to the second. 
In order to make close to 0 in the last line of Eq. (j^i — P2)tgr should be an odd integer 
multiple of n. Thus, the thinnest thickness, tg^^o is given by 
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Fig. 6. Transmittances of (a) the TM HCG and (b) the TE HCG as a function of wavevector, 
kx and angular frequency, ft). The green dashed lines designate 1550-nm wavelength, light 
line, and the 1st Brillouin zone boundary. The red dashed lines designate the guesses for 
the lowest guided modes below the light line. 


Substituting j 8 ^ values of the TE and TM waveguides given in Fig. 5 into Eq. ([3]) gives rgr,o 
values of 299 nm and 374 nm for the TE and TM HCGs, respectively. These values are similar 
as the exact values given in the caption of Fig. 1. The difference in the optimal thicknesses of 
TM and TE HCG reflectors is attributed to their difference in (j3i — ^ 2 )- 

This difference in (j3i — j 82 ) originates from the polarization-dependant boundary conditions. 
In Fig. 5, both TM mode 1 and TE mode 1 have more electric fleld conflned within the grating 
bars than the air, while both TM and TE modes 2 have more in the air. In the TM case, the 
continuity of Dx(=£Ex) is required at the interface of grating and air. As shown in Fig. 5(a), 
this results in a large enhancement of the electric fleld in the air region for the TM mode 1 
case, due to the large refractive index difference between the grating and air. However, in the 
TE mode 1 case, the continuity of Ey is required at the grating-air interface, which results in 
no enhancement as shown in Fig. 5(c). Since the propagation constant, j 8 ^ is proportional to 
the effective index of the mode, averaged with the electric fleld intensity as a weighing factor, 
the effective index of the TM mode 1 cannot be high: considerable electric fleld intensity is 
also distributed in the air while in other cases, the effective indices are determined solely by the 
Ailing factor. Therefore, (j3i — P 2 ) becomes smaller for the TM case. 

3.2. Grating period 

HCGs can be approximated as a slab waveguide with waves propagating along x direction. The 
transmittance plots in Fig. 6 show the dispersion of this waveguide, i.e., ft) vs. kx above the light 
line |[T9ll . Due to the periodicity, the waveguide dispersion is folded at the first Brillouin zone 
boundary, i.e., kx=7r/A. Since k^ ^ for the lowest waveguide mode (resonance condition 
along z direction), the kx value can be approximately evaluated in the unfolded dispersion curve 
by: 



In order to have a crossing with the ft) axis at (k^, ft))=(0, 2;r/1550 nm) in the folded dispersion 
curve as shown in Fig. 6 , Eq. (Hi should be satisfied at {kx^ ft))=(2;r/A, 2;r/1550 nm): 
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Fig. 7. Transmittance, (l—R) of the top reflector seen from the cavity layer as a function 
of the airgap thickness. A TE- and a TM-polarized plane wave at a wavelength of 1550 nm 
are incident from the cavity layer for the TE HCG and the TM HCG cases, respectively. 
The same grating parameters as in Eig. 1 are assumed. The refractive index of the cavity 
layer made of InP is 3.1661. 


( 5 ) 

Equation ^ shows that the larger period of the TE HCG can be attributed to its thinner thick¬ 
ness. Since the thinner thickness of the TE HCG originates from the boundary condition, one 
may say that the polarization dependent boundary condition leads to the polarization dependent 
differences in period as well as thickness observed in Eigs. 2(a) and 2(b). 

3.3. Evanescent tail 

The grating period determines this evanescent tail length. The propagation constant of the m-th 
order harmonic component in the air, is given by, 

kf = sJ{2K/Xf-{2nm/Af, (6) 

where A is the wavelength of an incident wave and A is the grating period. Since A < A, all 
higher order components than the zeroth decay evanescently. Equation ([b]) tells us that a smaller 
A gives a larger purely imaginary leading to a shorter evanescent tail. This explains why 
the TE HCG with a larger A has a longer evanescent tail. 

4. Discussion for HCG-based VCSEL designs 

As shown in Eig. 7(a), HCG reflectors require a low-refractive-index gap layer between a HCG 
layer and a cavity layer when incorporated into VCSELs. The reflectivity, R seen from the 
cavity layer is a composite reflectivity including the reflections from the HCG as well as the 
air gap/cavity interface. If the air gap thickness is close to the evanescent tail length, the HCG 
mode fleld may couple to the cavity layer, which reduces R. Eigure 7(b) compares the R values 
of the TM and TE HCGs of Eig. 1 as a function of the air gap thickness. The TM HCG reaches 
the saturated reflectivity value at an air gap thickness of 0.3 /im while the TE HCG does at 1.7 
/im. This is consistent with the observation that the TM HCG has a shorter penetration tail than 
the TE HCG. Note that the air gap thickness of 0.3 jim is 0.19 times of the wavelength, which 
is even thinner than the optical thickness of a single DBR layer. 

This thin gap thickness achievable for TM HCG allows TM HCG-based VCSELs to achieve 
a larger conflnement factor and a shorter effective cavity length than conventional DBR-based 
VCSELs. The larger conflnement factor reduces the threshold current, /th, increasing the mod¬ 
ulation speed through the dependence on stimulated emission current, (/ —This in- 







crease is more significant for low-enegy-consuming VCSELs for chip-level optical intercon¬ 
nects where the operating current, I is close to Ith. 1^ . Furthermore, the shorter effective cavity 
length is a key factor for achieving a higher modulation speed through a shorter photon life¬ 
time as well as for increasing the tuning range of HCG-based MEMS tunable VCSELs ID. For 
example, a TM HCG-based VCSEL can achieve a modulation speed higher than 100 Gb/s by 
using the findings of this paper ifTOll . 

The shorter period of TM HCGs is more adequate for wavefront engineering, such as beam 
focusing and tilting since it can make the phase front smoother. The reason is that the spatial 
variation of the reflection or transmission phase from a HCG is obtained by discrete phase 
variation from each period. 

As shown in Fig. 2(f), TE HCGs can be designed to have a long penetration depth, which 
is advantageous for achieving a large lateral coupling efficiency from the vertical cavity to an 
in-plane waveguide (91. The shorter penetration length of TM HCGs can allows for a short 
effective cavity length of HCG-based MEMS-tunable VCSELs, leading to a large tuning range 

M. 


5. Conclusion 

In conclusion, all of HCG parameters and some of HCG properties are different for the TM 
HCG and TE HCG broadband reflectors. Among HCG parameters, the grating thickness and 
grating period are quite different and these differences are found to originate from different 
electric field boundary conditions met for each type of HCG reflector. Among HCG properties, 
the reflectivity and stopband width are almost same for both types while the evanescent tail 
length and reflection phase are significantly different and the penetration depth is slightly dif¬ 
ferent for each type. Furthermore, it is found that the grating period and the grating bar width 
determine the evanescent tail length and the penetration depth, respectively. It should be noted 
that the polarization dependent HCG properties significantly influence the important properties 
of HCG-based VCSEL cavities including effective cavity length and quantum confinement fac¬ 
tor. Thus, they need to be considered in designing ultrahigh-speed VCSELs and MEMS-tunable 
VCSELs. 
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